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Abstract: We discuss the decay t" z/^7r^7r°7 in terms of a model with the 
correct low-energy structure and with the relevant resonance degrees of freedom. The 
nontrivial radiative dynamics becomes visible for large photon momenta only. We use 
the model to calculate electromagnetic corrections for the two-pion contribution to 
hadronic vacuum polarization extracted from photon-inclusive two-pion decays. The 
corrections are insensitive to the details of the model and depend on the pion form 
factor only. Putting all relevant isospin violating corrections together, we obtain a 
shift Aa^ = (—120 ±26 ±3) x 10~^^ to be applied to determinations of the anomalous 
magnetic moment of the muon from photon-inclusive r data. 
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1. Introduction 

The process r~ — > i/T-7r~7r°7 involves an interplay of strong, weak and electromagnetic 
effects. For small momentum transfer to the hadronic system as in the related radia- 
tive pionic beta decay vr"*" — Ueir'^'j , the theoretical determination of the associated 
decay amplitude is well under control. 

At low energies, the standard model of strong and electroweak interactions is de- 
scribed by an effective quantum field theory called chiral perturbation theory (CHPT) 
[0]. It treats the interactions of the experimentally observed hadronic degrees of free- 
dom at low energies by exploiting the constraints imposed by the symmetries of the 
standard model, in particular the spontaneously broken chiral symmetry. The ap- 
plication of CHPT is not restricted to purely strong interaction processes. With 
appropriate external source terms, also semileptonic and electromagnetic reactions 
involving leptons can be investigated. 

The applicability of CHPT ends once typical hadronic energies are approaching 
the mass of the p resonance. Nevertheless, the correct low-energy limit, unambigu- 
ously described by CHPT, puts constraints on the behaviour of amplitudes also in 
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the resonance region. The low-energy constants (LECs) of order in the chiral 
expansion are known to be dominated by the low-lying meson resonances 0. This 
correspondence facilitates the matching between the low-energy domain and the res- 
onance region. Applications along these lines comprise the pion form factor 0, 
radiative p decays and four-pion production 0. 

The radiative r decay is a further process where the resonance regime of the 
standard model can be tested. A second motivation for the present investigation lies 
in the relevance of the non-radiative decay Ut-tt^tt^ for a precise determination 

of hadronic vacuum polarization, which is essential for the calculation of both the 
energy dependence of the QED fine structure constant and of the anomalous magnetic 
moment of the muon. The reason is that a CVC relation between electromagnetic 
and weak form factors allows to relate the differential decay rate of r~ UrTc^Tc^ to 
the cross section a{e^e~ tt+tt") in the isospin limit. 

However, the present theoretical and experimental accuracy of hadronic vacuum 
polarization requires a precise determination of all isospin breaking effects due to 
rriu 7^ rrid and electromagnetism. The radiative corrections for r~ — > z/^7r~7r° are 
part of this task. To cancel the infrared divergences generated by the one-loop 
exchange of a virtual photon, the associated radiative decay mode must be included. 

In the low-energy limit, the necessary modifications of the CHPT framework in 
the presence of virtual photons || and virtual leptons have already been worked 
out. In this paper we want to cover also the resonance region relevant for the two-pion 
decay of the r lepton. Our present work is an extension of a previous letter |jlO[ where 
electromagnetic corrections and their impact on the determination of had already 



been discussed. However, the analysis of Ref. |T0[ carried a rather large uncertainty 
because of the incomplete matching between virtual and real photons that depends 
on the experimental setup. In the present paper, we discuss in detail the radiative 
corrections for the case where all radiative events are included in the two-pion data 
sample. This is a case of immediate practical interest because it corresponds to the 
analysis of the ALEPH experiment [pl]]^. 



The paper is organized as follows. In Sec. |^ we discuss the general structure of the 
— s> z/T-7r^7r°7 amplitude together with the dominant contributions in the resonance 
region. The decay rate is calculated in Sec. |^. Applying different cuts in the photon 
energy, we explore the energy region where deviations from bremsstrahlung can be 
expected. In Sec. ^we turn to the calculation of radiative corrections for the r decay 
into two pions relevant for hadronic vacuum polarization. The isospin violating 
corrections for are collected in Sec. ^. Finally, our conclusions are summarized in 
Sec. ^. The chiral resonance Lagrangian, kinematical details and loop functions are 
relegated to two appendices. 



^ We are grateful to Michel Davier and Andreas Hooker for information on the ALEPH procedure. 
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2. Amplitude for r i/^^ 7r% 

The matrix element for the decay 

has the general structure 

T = eGFV:,e'^iky{ FMq)Yil - 75)(m.+ f - ^hMP) (2-1) 
+ (V- VM?)7^(l-75)w(P)} . 

The first part of this matrix element describes bremsstrahlung off the initial r lepton 
with 

Fu = {Po - P-)J+{t)/2P -k, t = {p,+ poY . (2.2) 

The form factor governs the non-radiative decay. It is identical to the electro- 
magnetic form factor of the pion because we work in the isospin limit m„ = rrid- 

The second part of Eq. ( p.lD describes the vector and axial-vector components 
of the transition 

W-{P-q)^n-{p_)n\poMk) . 

The hadronic tensor V^i, contains bremsstrahlung off the vr" in the final state. Gauge 
invariance implies the Ward identities 

FV;, = (p_-po)./+W (2.3) 
k^-A^iy = . 

The tensor amplitudes V^v,^/xi/ contain four invariant amplitudes each |T^. For 
V^u we use a rather special form here that is convenient for our purposes. In A^^ we 
only include the two amplitudes that actually arise in our calculation: 

V,. = /+|(P - 'jf\f^(P- + * - Po). - M{P - lf\9^. (2.4) 

+ Mdf^uP- ■ k - p-^k^) + V2{gt,uPo ■ k - Po^k^) 

+ VsiP-f^Po ■ k - po^p^ ■ k)p^y + Vi{p^f,po ■ k - po^p^ ■ k){p^ +Po + k)^ 
Af,^ = iaie^^p^{po - p^yk" + ia2{P - q)u£f^paTk^p'Lpl ■ (2.5) 

These tensor amplitudes have the following properties: 

• Vfj,i, and A^i, satisfy the Ward identities ( P^) for any t>i, . . . , t>4, ai, 02. 
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Taking into account {P—qY = t+2{p-+pQ)-k, Low's theorem [13] is manifestly 
satisfied: 



V, 



p_ ■ k 



(2.6) 



+ 2 



p_ • k 

df+{t) ( p^^po ■ k 



dt 



P- ■ k 



-Pofi (P- -Po)u + 0{k) 



The invariant amplitudes vi,. . . ,V4 as well as 01,02 remain finite for A; — 0. 
Of course, initial state bremsstrahlung in Eqs. ( p.l| , p.2|) contributes also in the 
Low limit. 



The leading order in the low-energy expansion corresponds to 



/+ = 1 

Vi = V2 = Vs = V4, = 

Oi = 02 = . 



(2.7) 



Our aim is to construct a model for radiative r decay that can be trusted in all 
of phase space. A major constraint is the correct low-energy limit to O(p^). The 
local contribution of O(p^) contains information on which meson resonances must be 
taken into account for extending the chiral amplitude into the resonance region. On 
the other hand, except for the pion form factor f+{t), we are not going to include all 
one-loop contributions arising at 0{p^), let alone higher orders. Since we want to use 
the model for all energies accessible in r decay a complete one-loop calculation would 
not be of much help anyway. Part of the loop contributions will be contained in the 
energy-dependent widths of the resonance propagators and in the representation of 
the pion form factor of Ref. [0] that we shall use. Therefore, our approach uses 
the low-energy limit to identify the relevant degrees of freedom but the model goes 
clearly beyond the range of applicability of the low-energy expansion of QCD. Such 
a procedure works very well for the pion form factor itself p, ^ and it also seems to 
be successful in radiative p decays P and in four- pion production 0. 

In addition to the one-loop amplitude, V^^, depends to O(p^) on the LECs Lg, Liq 

m 

f+{t) = l + 2Ut/Fl 

vi = 4(2L9 + Lio)/F^ (2.8) 

f 3 = f 4 = . 
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Figure 1: Anomalous diagrams for the transition W it tt^^ contributing to the axial 
tensor amplitude A^^^. 

The constants Lg, Lio can be interpreted as low-energy remnants of p and ai exchange 
[0]. As usual, p exchange will turn out to dominate by far and the ai contribution is 
in practice undetectable in this decay. 

For completeness, we include also the (leading) contributions of 0{p^) to A^^, 
due to the chiral anomaly. The diagrams in Fig. |l| get contributions from the Wess- 



Zumino-Witten functional 14 



_ 1 _ 1 

~ 8^ ' ~ "47r2F2[(P-g)2-M2] ■ ^^-^^ 

By itself, the anomaly contribution cannot be expected to yield a realistic axial 
amplitude for all available energies. We take it rather as an indication of the expected 
order of magnitude of A^^. Anticipating the result of the numerical calculation, the 
anomalous amplitude is very much suppressed compared to the dominant p exchange 
contributions and turns out to be completely negligible in radiative r decay. 

The amplitude of O(p^) in ( |2.8| ) tells us that p and Oi exchange contribute to 
the vector amplitude at low energies. We now assume that these are the dominant 
degrees of freedom for all accessible energies. As is easy to understand from the 
decay rates involved, ai exchange will be much suppressed compared to p exchange 
except at very low energies. Nevertheless, we include the ai contributions to V^^. 

The resonance Lagrangian [0] and corresponding diagrams are displayed in App. |^. 
The Lagrangian gives also rise to double-p exchange diagrams. One of them is needed 
for ensuring gauge invariance beyond O(p^) accuracy but we include all of them for 
consistency. In terms of the couplings defined in App. the final result takes the 
form: 

n p 



+ - + t/^)D;:[iP- + kY] (2.10) 
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TT p 

+ - M'^-P- ■ k)D-,'[{p^ + m (2.11) 

TT ai 

-^ = TiLD-^'{ip-+kr] (2.12) 

-4 = -^-^D;\t)D;^[iP - qf] + ^,D;^{iP - g)1 • (2.13) 

TT TT p 

The resonance propagators are given by 

D^{t) = Ml-t-iMpVp{t) (2.14) 
Da,{t) = Ml^-t-iMaJaAt) 
in terms of the momentum dependent width [H] 



r,(0- 



9677^2 



AMl/tf''^e{t - AMI) + -(1 - AMl/tf^eit - AM] 



(2.15) 



For ra^(t) we use the parametrization of Ref. [|T5[. Because of the already mentioned 
suppression of ai exchange neither the form of (t) nor its on-shell value (we take 
Tai{M^J = 0.5 GeV) are of much importance in practice. 

To define the decay amplitude ( p.l| ) completely, we still have to specify the pion 
form factor f+{t). We use the results of Ref. where the chiral form factor of 0(j9^) 
was matched to the resonance region 

/+(t) = MlD;\t) exp \lH^^{t) + HKK{t)\ . (2.16) 

The subtracted loop function Hpp(t) is defined in App. 0. 

The pion form factor { '2.1t ) exhibits the correct low-energy behaviour to O(p^) 
by construction To check the low-energy limit of the vector amplitudes vi,V2 
(f3,f4 arise only at 0{p^)), one has to insert the resonance dominance relations p| 



^10 



AMI 



3. Decay rate and spectrum 

The differential rate for the decay r~ z/^7r"7r°7 is given by 

^ ' spins 
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BR(full) 


BR(brems) 


100 


8.4 ■ 10-4 


8.0 ■ 10-4 


300 


1.8 ■ 10-4 


1.5-10-4 


500 


3.8 ■ 10-5 


2.6 • 10-5 



Table 1: Branching ratios BR{t^ v^-ti^-k^^) for different cuts on the photon 

energy. The second column corresponds to the fuh amphtude of Sec. the third one to 
the complete bremsstrahlung approximation defined in the text. 



^^^^^ = 2£; 2k^2k:; 2F^ - ^ - ^- - - ^) • 

The decay rate is dominated by bremsstrahlung of soft photons. To be sensitive to 
the dynamics of the proper radiative transition, a substantial cut on photon energies 
will be necessary. 

We distinguish in the following between the full amplitude presented in the pre- 
vious section and what we shall call "complete bremsstrahlung", i.e. the amplitude 
with 

= ti2 = t>3 = t>4 = ai = 02 = . (3.2) 



Since we use the expression (p.l6| ) for the pion form factor in both cases, the correct 



Low limit ( |2.6|) is guaranteed in both scenarios. 

In Table |1|, we display the decay rate for three different lower cuts on the photon 
energy. In Fig. |^ we compare the t-distribution for > = 300 MeV (in the 
T rest frame) for the two options: the complete spectrum is shown as the full curve, 
the distribution corresponding to complete bremsstrahlung is given by the dashed 
curve. The interpretation of the spectrum is straightforward. The dominant peak is 
almost exclusively due to bremsstrahlung off the 7r~ . The secondary peak receives 
contributions from bremsstrahlung off the r lepton but for a sufficiently high cut on 
the photon energy the spectrum is also sensitive to the additional resonance exchange 
contributions in Vfj^^ due to the diagrams in Fig. |. 

There are at present no published results for the radiative decay under discus- 
sion. Table |I| and Fig. ^ suggest i?™"^ ~ 300 MeV as a reasonable value for the 
cutoff in photon energy to test the nontrivial dynamics of the radiative decay. For 
smaller values, rate and spectrum become more and more dominated by (complete) 
bremsstrahlung. For larger cutoffs the branching ratios become too small. 

We add a few comments on the theoretical status of the amplitude. Although 
our approach is well motivated from the low-energy structure predicted by QCD we 
cannot claim more than what it is worth: a realistic model for the intermediate-energy 
region that should be confronted with experiment. For instance, there is no guarantee 
that effects that only show up at 0{p^) and higher in the low-energy expansion are 
fully accounted for. The situation is different for complete bremsstrahlung. In this 
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Figure 2: Decay spectrum (arbitrary units) as function of t, the invariant mass squared 
of the two pions, for photon energies bigger than £'™™=300 MeV (in the r rest frame): full 
amplitude (full curve) vs. complete bremsstrahlung (dashed curve). 



case, the amplitude is completely determined by the pion form factor (|2.16|) which is 
known to describe the data very well in the region of interest 0. Therefore, whenever 
complete bremsstrahlung dominates (for E"™" < 100 MeV), rate and spectrum can 
be predicted in a model independent way. 



4. Radiative corrections for hadronic vacuum polarization 

The experimental precision in r decays has made very accurate determinations of 
hadronic vacuum polarization possible. In particular, a precision of 1 % has been 



achieved |jT6|, |lj] for the contribution of hadronic vacuum polarization to the anoma- 
lous magnetic moment of the muon. This accuracy makes a careful investigation of 
isospin violating and electromagnetic corrections mandatory. 



In a recent letter |[T^, we have calculated the electromagnetic one-loop correc- 



tions for the two-pion decay of the r lepton. Once the infrared divergences are 
removed by adding the decay rate due to soft photons, this loop correction is in 
principle well defined. However, a sizable uncertainty remains in general because 
the treatment of radiative events depends on the experimental setup. Without spec- 
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ifying those experimental conditions, we had to assign a rather big uncertainty to 
the radiative corrections. At best, the procedure of Ref. |jlO| may be apphed when 



only radiative events with photon energies below a sufficiently small energy cutoff 
are included. In this case, the leading Low approximation to the decay rate, as con- 
tained for instance in the widely used Monte Carlo program PHOTOS |T^, may be 



expected to provide a ffist approximation. But even in this case, the specific value 
of the cutoff in photon energy would still have to be implemented in the analysis of 



Ref. 10 



For the precision needed in the determination of hadronic vacuum polarization, 
a more careful analysis is mandatory. This is especially true for the case where 
all radiative events are included in the two-pion data sample as for instance in the 



ALEPH experiment [O]. When photons of all energies are considered the leading 
Low approximation for the radiative rate is of course no longer justified. 

For this purpose, we will use the model of Sec. ^ to calculate the radiative dif- 
ferential decay rate. To make the cancellation of infrared divergences transparent, 
we divide the calculation into two parts. The infrared divergent decay rate in lead- 
ing Low approximation will be calculated analytically. When added to the equally 
infrared divergent electromagnetic one-loop correction [|r^ for the non-radiative de- 
cay, the infrared divergences cancel and the limit of vanishing photon mass can be 
taken. The second part of the radiative rate, containing in particular the subleading 
Low terms of 0{k~^), will be calculated numerically. This second part is infrared 
convergent and the photon mass can be set to zero from the start. 

In the notation of our letter [^, we set out to calculate the electromagnetic 
correction function G'em(^) defined by 

%*^ = ^^i#^(i-«4V«P(i-ij«i + ^)IAWPt?EMW . (4.1) 

The decay rate here is the fully inclusive one, i.e., with photons of all possible energies 
included. S'ew is a short- distance electroweak correction factor [|19|. We also recall 
that we work in the isospin limit in this section. 

In a first step, we consider only the leading Low approximation of 0{k~'^) to the 
differential decay rate for the radiative process: 

dT='^^^^\U{t)?D{t,u)x (4.2) 



327r'^mT- 

2P ■ p_ ml Ml 



{P ■ k - My2){p. ■ k + MV2) {P-k-MV2y (p_-fc + M2/2)^ 



d 



LIPS • 



The kinematics of the radiative decay and the function D{t,u) are given in App. |B|. 
The Lorentz invariant phase space element dups is defined in Eq. ( |3.1| ). Here and in 
the following, we neglect terms that will not contribute to the final result for — ^ 0. 
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Integration over neutrino and photon momenta leads to the three-fold differential 

rate 



dT 




'^\U{t)\'D{t,u)x 



(4.3) 



{2P ■ p_Iii(t, u, x) — mll2o{t, u, x) — M^/o2(t, u, x)} dt du dx 



with 




n 



(4.4) 



In the next step we perform the integration over x, the invariant mass squared of 
photon and neutrino. Here we must distinguish between two different regions in the 
t — u plane. For (t, u) in the Dalitz plot of the non-radiative decay, the photon mass 
must be kept nonzero and the lower limit for the x-integration is in fact M^. For 
{t, u) that cannot be accessed in the non-radiative decay the lower limit is given by 
X-{t,u) in (p.5| ). The resulting contribution to the rate is infrared finite. The latter 
case can only occur for t < m^MT^/ {nij. — M^) = 0.27 GeV^. The corresponding 
contribution to the spectrum is very small. We will therefore not display the explicit 
formula for this part of the spectrum here but it will be included in the correction 
function Gem(^)- The upper limit of the x-integration is always given by x+(t,-u) in 



The double differential decay rate in the leading Low approximation, for (t, u) 
in the non-radiative Dalitz plot and integrated over all photon momenta, takes the 
form 



Neglecting as always terms that vanish for 0, the functions Jmn are given by 



dT 




|L|/4t)|2/)(t,^)x 



(4.5) 



{Jii(t,M,M^) + J2o(t,M,M^) + Jo2it,u,M^)} dt du . 



Ju{t,u,M^) = log( 




+ i {Li^il/Y^) - Li^iYi) + log2(-l/r2)/4 - log2(-l/Fi)/4} 




Jo2(t,u,M^) 



(4.6) 



The various quantities in these functions are defined in App. 0. 

The remaining part of the radiative decay rate, containing all terms except the 
leading Low terms, is calculated numerically. Putting everything together, the double 



differential rate for the photon-inclusive two-pion decay is given by 



dt du QAir'^irfi 



\U{t)\^D{t,u) {l + 2f^X{u,M,)+g,Ut,u,M,)] . (4.7) 



The electromagnetic loop amplitude /jq™(u,M^) was calculated in Ref. [jTO| and is 
reproduced in App. The radiative part of the rate is determined by the function 

gri,d{t,u,M^) with 

grsdit, U, M^) = 5'brems(t, U, M^) + fi'rest(^, u) (4.8) 

a 

gbrems{t,U,M^) = -[Jii{t , U, M^) + J2o{t,U,M^) + Jo2{t,U,M^)] . (4.9) 

The (numerically calculated) function grest{t,u) accounts for the infrared finite re- 
mainder of the rate, i.e., all terms except the ones contained in (7brems(^, ^, ^7)- As 
it must be, the sum + ghrems is independent of the photon mass and therefore 

infrared finite as well. 



To obtain the correction function G'em(^); we still have to integrate over u [|T0 

(t) 

duD{t, u 



duD{t,u) A{t,u) 
GEuit) = , (4.10) 



A{t,u) = l + 2f^^(u, M^) + ^?brems(t, M^) + g,,,,{t, u) . (4.11) 



with 

^ '^^ "'^ ^ ' /loop I "'^''7/' "I" i/bremsV''; '^5 -"'■'7 

The limits of integration Umin,max(^) can also be found in App. 0. 

Since only the sum -f^fbrems is physically meaningful two distinct pieces can 
be distinguished in G'em(^)- Setting g^^^i to zero, we obtain a function G^^l^{t) that 
describes radiative corrections in the leading Low approximation. In Fig. |^ we plot 
both the complete result Gem(^) and the (leading) bremsstrahlung approximation 

The obvious conclusion is that the non-leading terms in the Low expansion can- 
not be neglected when considering the photon-inclusive spectrum and rate. One 
should keep in mind that the function Gem(^) multiplies the non-radiative spectrum 
in ( [4. ID that is strongly peaked around t = ^0.6 GeV^. 

Where does the difference between the two curves in Fig. ^ mainly come from? 
The answer is that it comes almost exclusively from the subleading Low terms in 
(12.61). In other words, the amplitudes f 1, . . . , ^4, ai, 02 have practically no impact on 
G'em(^)- In fact, G'em(^) for 

vi = V2 = V3 = V4 = ai = a2 = 
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1.04 




0.5 1 1.5 2 2.5 3 

t [GeV^] 

Figure 3: Correction function GemIO defined in (|4.l| ) and ( [4.10D . The full curve corre- 
sponds to the complete radiative amplitude, the dashed curve describes the leading Low 
approximation . 



could not be distinguished from the full result in Fig. |^. This conclusion corresponds 
to the findings in Sec. |^ that the specific features of the model for the radiative 
transition only show up for a sufficiently large lower cut on photon energies. 

In the present case where photons of all energies are included, complete brems- 
strahlung is an excellent approximation. However, as Fig. |^ documents very clearly, 
the leading Low approximation (0(A;~^) in rate) corresponding to the function G^l^{t) 
is not a valid approximation. The complete bremsstrahlung approximation {vi = 
V2 = = = ai = a2 = in (|2.4|), (p.5[)) is uniquely determined by the pion form 
factor f+{t). The correction function GemI^) is therefore effectively independent of 
the details of the model of Sec. for the radiative transition. 

Finally, we emphasize once again that the analysis in this section applies to the 
photon-inclusive case only. Moreover, even the leading Low approximation G^l^{t) 
(dashed curve in Fig. |]) cannot be compared directly with the corresponding function 
in Ref. |jlO[ that contained only the loop contribution (with an ad hoc prescription for 
the infrared divergent part). In contrast, the curves in Fig. prefer to the complete 
electromagnetic corrections, i.e., to the (infrared finite) sum of loop and radiative 
contributions with photons of all energies included. 
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5. Isospin violating corrections for a 



p 



The leading hadronic contribution to the anomalous magnetic moment of the muon 
O'lj. = (A'm — 2)/2 is due to hadronic vacuum polarization pO[ , 



vacpol 



47]-3 



^hadrons 



(5.1) 



denotes 



where K{t) is a smooth kernel concentrated at low energies and c"e+e-^hadrons 
the bare hadronic cross section with QED corrections removed. The low-energy 
structure of hadronic vacuum polarization is especially important. In fact, about 70 



% of a^'^P"' is due to the two-pion intermediate state for 4M^ < t < 0.8 GeV' 
Ref. 



[see, 

e.g., Jrtei. |i21|). A substantial improvement in the accuracy is possible if one includes 
data from hadronic r decays. A precise link between hadronic spectral functions from 
r decays and from the e^e~ hadronic cross section requires the calculation of radiative 
corrections as well as the inclusion of other isospin breaking effects (both kinematical 
and dynamical in origin). For the two-pion final state, we have parametrized |]10| 
the relation between the bare e^e~ cross section and the observed differential r decay 
rate as follows ^: 



Krit) dt 



X 



Rmit) 
Sew 



(5.2) 



where 



Kr{t) 



iKdpm? 



vr a 



3t 



3847r3 

2 



1 + 2- 



(5.3) 



and the isospin breaking correction is in 



Rm{t) 



1 Pi 
Gem(^) Pi 



Fv{t) 



hit) 



(5.4) 



The factor Sesn 



1.0194 (taken at the scale m^) takes into account the dominant 

Rmit) involves the long-distance QED 



short-distance electroweak corrections 19 



correction G'em(^) (discussed in detail in the previous section), the phase space cor- 
rection factor (t)//3j|()jr- (^) (important in the low-t region), and the ratio of 
electromagnetic {Fy) over weak (/+) pion form factors. 

In this section we discuss the isospin breaking affecting the ratio of form factors. 
Since our purpose is to relate the observed weak decay to the bare two-pion cross 



^The bare e+e~ cross section directly provides the two-pion contribution to the muon g — 2 at 
O(a^). We follow here the tacit convention that only explicit photonic corrections are considered 
in the counting of powers of a. In other words, electromagnetic contributions to the charged pion 
mass or to p-Lo mixing are included in the lowest-order contribution of 0{a^). 
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section, we need to include effects of order (m„ — rud) in both Fy and /+, while we 
need explicit photonic corrections of order a only in /+. Keeping track of all relevant 
sources of isospin breaking at first order, Eq. ( |2.16| ) generalizes to 



Fy{t) = M\D-}{t) 



exp [2H^+^-{t) + HK+K-{t) 



9 



pu) 



t 



U{t) = M'D.\t) exp 2H^^At) + HK+Mt) + f^T.i + ■■■ 



3M2M2 

pelm 



t - 



(5.5) 
(5.6) 



In Eq. ( |5.5|) the main new ingredient compared to Ref. |]T0[ is the inclusion of the u 
contribution to the electromagnetic form factor, following the notation of Ref. [p^]. 



This contribution was not included in our previous work ||T0[ because it is formally of 
higher order in the low-energy expansion. Here we include p-u mixing together with 
all other numerically relevant isospin violating contributions. The mixing parameter 
can be extracted from the decay width for uj — > vr+vr" leading to 9p^ = (—3.9 ± 
0.3) ■ 10~^ GeV^ in the framework of Ref. p2[. A more direct determination comes 
from the recent Novosibirsk data for e~^e~ tt+tt^ |^^. Translating their fit for 
the mixing parameter into our notation gives 6pi^ = (—3.3 ± 0.5) ■ 10^'^ GeV^ (the 
same value has been obtained in an independent fit of the CMD-2 data p4|). For 
the numerical discussion we take 6 



piU 



-3.5 ± 0.7) ■ 10 ^ GeV^. For the remaining 



parameters we use M, 
Eq. 



0.783 GeV and F,, = 0.00844 GeV. 



]6D reproduces our previous result ||T0| where represents the structure 



dependent electromagnetic corrections to the low-t part of the spectrum: 

.2 n/r2 ^2 



^elm 
J local 



a 
Att L 



3 1 m; 



log^ + 21og^ 



-(4vr) 



1 



-2KIM + -X, + -x^M 



(5.7) 



Ki2, Xi and Xq are LECs appearing in the effective Lagrangian of CHPT. Numerical 
estimates for the chiral couplings are reported in Ref. 



IG]. The dots in Eq. (E 



represent unknown structure dependent corrections to the form factor, away from 
the threshold region. At the moment there is no model independent method for 
including such terms. However, our calculation correctly identifies the large UV logs 
(in Sew) and the IR logs (chiral logs in f[^^i and /lo^), and therefore the neglected 
terms are expected to be numerically smaller. 

In the numerical analysis we also account for isospin violating effects due to 
mass and width differences between the charged and the neutral p mesons. The p^- 

We follow Ref. ^ 



mass difference is unknown but expected to be small [^, [TT 



16 



and assume Mp+ — MpO = (Oil) MeV. In addition to the two-pion contributions, 
there is also a possible width difference due to radiative decays, dominated by the 
nnj modes. We include the measured TTir'-f decay width of the p° of approximately 1.5 
MeV ||26| and assume a radiative decay width difference T[p^ tt^ti^^) — F(p+ 
TT+TT^) = (0.45 ± 0.45) MeV fS^, |, [T§. 
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Figure 4: The ratio for = (— 3.5±0.7)-10^'^ GeV^. The full curve corresponds 

to the mean value. For this plot we have assumed Mp+ = M^o = 775 MeV and we have 
neglected the radiative decay widths. 



We can summarize our results graphically. In Fig. ^ we plot the ratio 
showing the sensitivity to dp^. In Fig. ^ we plot the full isospin violating correction 
function -Rib(^) containing the complete radiative correction GEuit), displaying also 
the small uncertainty due to the electromagnetic LECs. 

In order to investigate the impact of RiB(t) on a™^P°\ we have evaluated for the 
two-pion contribution to the difference between the correct expression for a^^ 
(full RHS of Eq. ( p.2| )) and the naive CVC expression (the term in square brackets 
on the RHS of (|5.2| )), up to a scale tmax: 



^^vacpol 



47r3 



4M2 



dtK{t) 



r[7] 



Kr(t) dt 



Rm{t) 



EW 



(5. 



In Table ||, we have analysed separately the effects of S'ew and of the three ingredients 
of Rib'- kinematics (KIN), radiative correction (EM) and form factor ratio (FF). We 
find that the difference Aa^'^P"' is quite insensitive to tmax as long as tmax > 1 GeV^ 
(cf. Table 0). We can summarize our results as follows: 



The short-distance correction S'ew 



1.0194 induces Aa™' 



-97-10 



-11 
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Figure 5: Full correction factor Rmit) with the complete radiative correction Gem(^) 
(solid curve). The (almost indistinguishable) dashed curves correspond to variation of the 
chiral couplings within their error bars |jl^. In this plot we use 9puj = —3.5 • 10~^ GeV^, p 
masses and widths as in Fig. ^. 





Sew 


KIN 


EM 


FF 




1 


- 95 


- 75 


- 11 


61 ± 26 ± 3 


- 119 


2 


- 97 


- 75 


- 10 


61 ± 26 ± 3 


- 120 


3 


- 97 


- 75 


- 10 


61 ± 26 ± 3 


- 120 



Table 2: Contributions to Aa™'^^" from various sources of isospin violation (in units of 
10^^^) for different values of tmax (in units of GeV^) as defined in Eq. ( ^.81) . 

• The kinematical effect induces a shift Aa™'^^"' = —75 • 10~^^. 

• The complete radiative corrections produce a shift Aa™^P°^ = —10 • 10~^^. We 
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emphasize once more that this correction apphes only for the case of a photon- 
inclusive two-pion data set as provided for instance by the ALEPH experiment 
|TT| . If we had used the leading Low correction factor (^^^(t), the shift would 



be Aar^P°i = +16 ■ 10"^^ instead. 



The ratio of form factors produces a positive shift of approximately 61 ■ 10"^^. 
Roughly 60 % of this effect is due to p-u mixing, the remainder comes from 
the p^-p° width difference (pion mass difference and radiative decay modes). 
We add the errors due to the previously discussed uncertainties of p-u mixing 
(8), p~^-p~ mass difference (20) and radiative decay widths (14) in quadrature 
(the numbers in brackets denote the corresponding shifts in units of 10^^^), 
but keep the small additional uncertainty in the electromagnetic LECs |jT^ (3) 
separate. The final shift due to the form factor ratio is therefore /\aj^^^°^ = 
+ (61 ±26 ±3) ■ 10-". 



Our results are based on the form factor representations of Eqs. ( |5.5| ) and ( |5.6| ). 



Use of the experimentally measured /+, as well as a different treatment of the form 
factor ratio [|l6l, produce slightly different values for the individual contributions to 
^^vacpoi^ rjij^g overall picture emerging from the two treatments is however quite 
consistent: there is a strong cancellation of all effects included in Ri^{t). The results 
in Table ^ put this statement on a quantitative basis. 

Collecting all isospin violating contributions, we obtain finally 

Aa™' = (-120 ± 26 ± 3) ■ 10"" . (5.9) 

We have not assigned an error for the short-distance correction factor Sew because 
this is already included in the error for the electromagnetic LECs. Sew is in fact 
closely related to the LEC Xg(Mp) after scaling it down from rrir to Mp |]10|, ^ 



Altogether, we are confident that with the photon-inclusive two-pion decay of the r 
as input the isospin breaking shift ( |5.9|) can be given to an accuracy of about 20%. 
The accuracy could be improved considerably by better determinations of the p mass 
and radiative decay width differences. 



6. Conclusions 

Guided by the low-energy limit dictated by CHPT, we have constructed a meson res- 
onance model which allows the investigation of the radiative r decay t~ — > Urn'ir^'j 
for the whole range of momentum transfer to the hadronic system. We have shown 
that for a photon energy cut of approximately 300 MeV in the r rest frame, the 
nontrivial part of the amplitude can be discriminated from the pure bremsstrahlung 
contribution in the decay spectrum. 

As a second application of our model, we have considered the radiative correc- 
tions for the process z/^tt^tt'^. This includes also real photon emission via 
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the associated radiative decay. Such an investigation is relevant once experimental 
data from hadronic r decays are used for the determination of the hadronic vacuum 
polarization, which requires a reliable estimate of isospin violating effects. 

The contribution from real photon emission to the decay rate depends on the 
experimental setup. If photons of all energies are included (as in the ALEPH data 
PT|), the leading Low approximation is definitely not sufficient any more. On the 
other hand, complete bremsstrahlung turns out to be an excellent approximation in 
this case. Compared to the naive CVC relation, we find a total isospin violating 
shift of the anomalous magnetic moment of the muon of Aa^'^P"^ = (— 120 ± 26 ± 
3) ■ 10~^^. The first error is due to the combined uncertainties of p-uj mixing, p^-p^ 
mass and (radiative decay) width differences. The comparatively small second error 
is based on available estimates of the low-energy constants. 
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A. Resonance contributions 

The resonance exchange contributions considered in Sec. |^ are derived from the 
standard chiral resonance Lagrangian [0] for vector and axial-vector mesons (scalar 
exchange does not contribute to our decay): 

C[V{l--),A{l++)]= C^JyV,A] (A.l) 

The Lagrangian is formulated in terms of SU (3) tensor fields V^u, (not to be 
confused with the tensor amplitudes of Sec. We refer to Ref. for definitions 
and further details, in particular for the explicit form of the gauge and chiral invariant 
kinetic terms that contain the couplings bilinear in the resonance fields. 

In addition to Mp = 0.775 GeV, M^, = 1.23 GeV and P^^ = 0.5 GeV we use for 
the resonance couplings the theoretically favoured values 



V2F^ = 0.13 GeV, Gy = F^/V2 = 0.065 GeV, Fa = = 0.0924 GeV. 



These values for Fy, Gy are used in the energy dependent width (|2.14| ) leading to 
Tp = Tp{M^) = 0.147 GeV. 

The vector amplitudes f i , . . . , are derived from the diagrams in Fig. 13. 
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Figure 6: p and ai exchange diagrams for the model of Sec. ^. In the first diagram, a 
photon is to be appended wherever possible. 

B. Kinematics and loop corrections 
B.l Kinematics 

After integration over neutrino and photon 4-momenta, the spin-averaged decay dis- 
tribution for 

depends on three invariants^. We use 

t = {Pi+P2)\ u={P-pi)\ x={q + kf (B.l) 

with = m^, pf = {i = 1,2). In order to calculate differential rates and 
spectra, we need the physical region V in the form of a normal domain: 

^ {^min ^ t ^max; Wmin(^) ^ ^ ^ ^max(^)) •^min (t,M) < X < Xmax(t,M)} , (B.2) 

^In the body of the text we use the more suggestive notation Pi ^ P-, P2 ^ Po- 
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where 



imin = (Ml + M2) 



(B.3) 



In order to give the exphcit form of Umin/max(^) and Xrain/max(t,u) we introduce the 
auxihary functions 



2t 



2t{ml + Mf) - (t + - M^) {ml + t~ x) 



± ^A(t,M2,M|)A(m2,t,x) 



(B.4) 



a;±(t,M) 



2M/ 



2 M/ (m^ + t) - (t + - M|) (m^ + M/ - u) 



± JA(t,M2,M|)A(w,M2,m2) 



(B.5) 



with A(a;, y,z) = x'^ + y"^ + — 2 {xy + xz + yz). In the following, let us call TZ^^^ 
the region in the t-u plane accessible in the non-radiative three-body decay, and 7^^^ 
the region accessible in the radiative decay. TZ^^^ is given by 



Mmin(t) = M-(t,0) 



for tmin ^ t ^ ^max 
for ^min ^ t ^ ^max 



(B.6) 
(B.7) 



while TZ^^ corresponds to 



Mr 



for tniin <t <U 



(B.8) 



,M+(t,0) 



for U<t < trai,x 
for tmin <t<t^ 



(B.9) 



with 



rririmrMi + M| - 



- Ml 

Finally, for given (t, m) the limits for x are 



(B.IO) 



(B.ll) 



M2 for (t, u) e TZ^^^ 

x-{t,u) for {t,u) en^^ 



The differential decay rate in Eq. (^.5|) involves the function 



2 

D{t, u) = ^{ml -t) + 2MlMl - 2u {ml - t + + Ml) + 2u^ . (B.13) 
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Moreover, the following combinations of invariants appear in Eq. 



1 - 2a ± v/(l - 2a)2 - (1 - /52) 



with 



a 



7 



{ml - t){ml + Ml - t - u) X{u,Mf,ml) 

Ml + ml-u 
^\{u,Mlml) 

-M^ml + Ml (m^ - t){Ml - u) - tu{-ml + t + u) 
+M^(-mt + tu + mlt + mlu) 



B.2 Loop functions 

The loop function Hpgit,^) is given by 



—X{t,Ml,M')r^{t) 



1 

12 



_12t 



A 



PQ 



' 18(47r)2 

(Apifi) - Agifi)) - 2(Ap(/i) + Agif,)) 



where 



EpQ = Ml + Mj 



A 



PQ 



Ap(/i) 



M|. 

'(47r)2 
1 



log 



Ml 



327r2 



2 + — ^ log ^ 



M2 



Mi -Mi 



Sp, M| 
^M2 



AV2(t,M2,M2) 



log 



ApQ 

(t + AV2(t,M2,M2))2-A 2.Q' 
(t-AV2(t,M2,M2))2_A 



PQ 



(B.14) 



,(B.15) 



(B.16) 



The one-loop virtual photon correction, contributing to the differential spectrum 
, is given by ^ 



a 
Air L 



[u - Ml) A{u) + {u-Ml- ml) B{u 

M^^m-r' 



+ 2{M^ + ml- u) C{u, M^) + 2 log ■ 



M2 



(B.17) 



^The last term in fi^^piu,M^) was missing from Eq. (16) of Ref. [ p!o| . This typo did not affect 
the numerical analysis of Ref. pOl . 
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In terms of the variables 

,2 



''^ = ^^ = 1 + - ' = Y7r^^^ ~ Vl/? - 4rr) , (B.18) 



and of the dilogarithm 

Lhiz)^- [ ^log{l -zt) , 
Jo I 

the functions contributing to /io^('W, My) are given by 



(B.19) 



A{u) = - 
u 



Bin) 
C{u, M^) = 



u 



2 ^ 1 - 
- log Tj. H ^ log Xr 



2 

1 x. 



r,- I — xi 



(B.20) 
(B.21) 



nirM^ 1 — xi 



1 TT^ 1 

-- log^ + 2 logX^ log(l ~ ^r) - Y + g log^ 



T M 



2 H 



+ Li2{x^) + Li2(l + 1/^2(1 - Xr-s/r^) - log log 



. (B.22) 



- 22 - 



References 

[1] J. Gasser and H. Leutwyler, Ann. Phys. 158 (1984) 142; Nucl. Phys. B 250 (1985) 
465. 

[2] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311. 

[3] F. Guerrero and A. Pich, Phys. Lett. B 412 (1997) 382. 

[4] D. Gomez Dumm, A. Pich and J. Portoles, Phys. Rev. D 62 (2000) 054014; 
J.A. Oiler, E. Oset and J.E. Palomar, Phys. Rev. D 63 (2001) 114009. 

[5] K. Huber and H. Neufeld, Phys. Lett. B 357 (1995) 221. 

[6] G. Ecker and R. Unterdorfer, Eur. Phys. J. C 24 (2002) 535. 

[7] R. Urech, Nucl. Phys. B 433 (1995) 234. 

[8] H. Neufeld and H. Rupertsberger, Z. Phys. C 68 (1995) 91; Z. Phys. C 71 (1996) 
131. 

[9] M. Knecht, H. Neufeld, H. Rupertsberger and P. Talavera, Eur. Phys. J. C 12 
(2000) 469. 

[10] V. Cirighano, G. Ecker and H. Neufeld, Phys. Lett. B 513 (2001) 361. 

[11] R. Barate et al. (ALEPH), Z. Phys. C 76 (1997) 15; Eur. Phys. J. C 4 (1998) 409. 

[12] J. Bijnens, G. Ecker and J. Gasser, Nucl. Phys. B 396 (1993) 81. 

[13] F.E. Low, Phys. Rev. 110 (1958) 974. 

[14] J. Wess and B. Zumino, Phys. Lett. 37 B (1971) 95; 
E. Witten, Nucl. Phys. B 223 (1983) 422. 

[15] J.H. Kiihn and A. Santamaria, Z. Phys. C 48 (1990) 445. 

[16] R. Alemany, M. Davier and A. Hocker, Eur. Phys. J. C 2 (1998) 123; 
M. Davier and A. Hocker, Phys. Lett. B 435 (1998) 427; 

A. Hocker, Talk given at "XXXVI Rencontres de Moriond", Les Arcs, March 2001, 
hep-ph/0111243 . 

[17] M. Davier and A. Hocker, analysis in preparation. 

[18] E. Barberio and Z. Wijs, Comp. Phys. Comm. 79 (1994) 291. 

[19] W.J. Marciano and A. Sirlin, Phys. Rev. Lett. 61 (1988) 1815; ibid. 71 (1993) 3629. 

[20] M. Gourdin and E. de Rafael, Nucl. Phys. B 10 (1969) 667. 

[21] S. Narison, Phys. Lett. B 513 (2001) 53. 



- 23 - 



[22] R. Urech, Phys. Lett. B 355 (1995) 308. 

[23] R.R. Akhmetshin et al. (CMD-2), Phys. Lett. B 527 (2002) 161. 

[24] J. Portoles, Talk given at "QCD 02", Montpellier, July 2002, to appear in the 
Proceedings. 

[25] J. Bijnens and P. Gosdzinsky, Phys. Lett. B 388 (1996) 203. 

[26] Review of Particle Physics (Particle Data Group), Eur. Phys. J. C 15 (2000) 1. 

[27] P. Singer, Phys. Rev. 130 (1963) 2441; Err. 161 (1967) 1694. 

[28] V. Cirigliano et al., Eur. Phys. J. C 23 (2002) 121. 

[29] G. Ecker et al., Phys. Lett. B 223 (1989) 425. 



- 24 - 



